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Abstract. Based on first-principles electron structure calculations and 
employing the frozen-magnon approximation we study the exchange interactions 
in a series of transition-metal binary alloys crystallizing in the zinc-blende 
structure and calculate the Curie temperature within both the mean-field 
approximation (MFA) and random-phase approximation (RPA). We study two Cr 
compounds, CrAs and CrSe, and four Mn compounds: MnSi, MnGe, MnAs and 
MnC. MnC, MnSi and MnGe are isovalent to CrAs and MnAs is isoelectronic with 
CrSe. Ferromagnetism is particular stable for CrAs, MnSi and MnGe: All three 
compounds show Curie temperatures around 1000 K. On the other hand, CrSe and 
MnAs show a tendency to antiferromagnetism when compressing the lattice. In 
MnC the half-metallic gap is located in the majority-spin channel contrary to the 
other five compounds. The large half-metallic gaps, very high Curie temperatures, 
the stability of the ferromagnetism with respect to the variation of the lattice 
parameter and a coherent growth on semiconductors make MnSi and CrAs most 
promising candidates for the use in spintronics devises. 



PACS numbers: 75.47.Np, 75.50.Cc, 75.30.Et 

1. Introduction 

The emerging field of spintronics [21 , also known as magneto- or spin-electronics, 
has provided a strong motivation towards the fabrication of novel magnetic materials. 
A strong interest has been focused on the so-called half-metallic ferromagnets, i.e. 
metallic compounds for which either the majority or minority spin-channel presents 
a gap at the Fermi level E F Such behavior is detected in various perovskite 
structures 0|, some oxides like CrC-2 [5] and magnetite (Fe304) [Sj, several Heusler 
alloys or some diluted magnetic semiconductors (DMS) [TulITT] . 

In 2000, Akinaga and collaborators succeeded to grow a thin film of the binary 
alloy CrAs on top of GaAs ^2- It was found that CrAs was actually continuing 
the zinc-blende (zb) structure of GaAs crystallizing in a mctastable phase. The 
most astonishing property was the total spin magnetic moment of 3 fiB, which is 
consistent with the half-metallic character predicted by relevant electronic structure 
calculations published in the same paper and reproduced with larger accuracy in latter 
publications ^3 • The thin films of CrAs have also the great advantage of a high 
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Curie temperature Tq, exceeding the 400 K (limit of the magnetometer used in the 
experiment) JS]. This discovery intrigued the interest of experimentalists and several 
attempts were made to grow similar metallic compounds as thin films. Zhao et al. have 
grown zb CrSb films by molecular beam epitaxy and found a Tq similar to the CrAs 
case Moreover, the growth of small zb MnAs islands on GaAs substrates ^7] 

has been achieved and CrAs/GaAs multilayers have been fabricated ^B] showing that 
coherent heterostructures are possible. These half-metallic materials combine high 
values of Tc with the coherent growth on top of semiconductors and as the know- 
how on the growth increases, they can be considered as very promising candidates for 
spintronics applications. 

Several theoretical studies have been devoted to the zinc-blende compounds 
of transition elements with group-V and VI elements (respectively, pnictides and 
chalcogenides). The calculations of the total energy and structural optimization 
show, in agreement with experiment, that the ground state crystal structure of these 
compounds is not of the zinc-blende type. In most cases the equilibrium crystal 
structure is of the NiAs type ^5]. This structure does not lead to half-metallicity 
|2"H1 PHI |22 12H1 E2| ■ The zinc-blende structure can only be obtained by the epitaxial 
growth of few layers on semiconductors. Galanakis and Mavropoulos studied in 
detail the magnetic properties and the position of the Fermi level in zb V-, Cr- and 
Mn-pnictides and chalcogenides as a function of the lattice constant [23 ■ Several 
authors have confirmed these results |26l 1271 1281 1291 1301 1311 13"2"] . Some of the binary 
semiconductors crystallize in the wurtzite structure and a number of studies is devoted 
to such systems j^Hl E]- Of special interest is the case of half- metallic zb MnC 
compound whose properties differ from the properties of other Mn compounds |35| . 
Also the cases of unsupported clusters [201 and nanocrystallites jjjjZI have been studied. 

An important aspect of the studies is the stability of the half-metallic gap. It 
has been shown that the tetragonal deformation of the zb lattice does not destroy 
the gap OH] and for some systems it even makes the zb structure more stable 
1 101 141j . Contrary to other half- metallic systems, in the case of the zb transition 
metal compounds it is possible to retain the half-metallicity both at the surfaces of 
these materials 123 021 and in the case of multilayers with binary semiconductors 
US EH ESH1I3 El EH]. The V and ^ r atoms have the ability to increase considerably 
their spin moment as a result of decreasing number of neighbors keeping in this way 
the half-metallicity of the system. 

Contrary to Cr and V, the Mn spin moment is saturated in these compounds that 
leads to the loss of the half-metallicity at the surface [2H1- The spin-orbit coupling 
has very weak effect on the spin-polarization of the electron states at the Fermi level 
for the sp atoms up to 5th period of the Mendeleev's table In the case of MnBi, 
where Bi is a heavy element from the 6th period, the spin-orbit coupling leads to 
the decrease of the spin-polarization from 100% to 77% ^Hj- Recently, Chioncel and 
collaborators studied the influence of the correlation effects on the electron structure 
of CrAs [SUj. They found that the spin-magnon interaction leads to the appearance 
of nonquasiparticle states in the spin-minority channel. The states are shown to lie 
above the Fermi level and to be sensitive to the lattice constant. 

In contrast to the gap properties, the exchange interactions in these compounds 
attracted much less attention. Shirai has shown that in the case of pnictides (VAs, 
CrAs and MnAs) the ferromagnetic state is energetically preferable compared to the 
antiferromagnetic state |13J . Only in the case of zb FeAs an opposite trend was found 
[T3] . Sakuma demonstrated that the Mn compounds isoelectronic to CrAs (MnSi, 
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MnGe and MnSn) are half-metallic ferromagnets |51|, In these compounds, similar to 
CrAs, the long-range ferromagnetism is stabilized by the short-range direct Mn-Mn 
interaction and an indirect interaction through the sp atoms. The Curie temperature 
in both CrAs and the isoelectronic MnZ compounds was found to exceed 1000 K 
|51|. Kiibler calculated Tq for VAs, CrAs and MnAs and found the largest value of 
1041 K in the case of CrAs compared with 764 K in VAs and 671 K in MnAs 
Taking MnAs as an example, it was shown that the loss of the half-metallicity in a 
compressed lattice leads to a strong decrease of the Curie temperature |S2j- Sanyal 
and collaborators have calculated Tq for V-, Cr- and MnAs compounds assuming the 
GaAs and InAs lattice constants |28|. Their mean field approximation (MFA) results 
agree with the calculations of other authors while the values of Tq obtained with 
Monte Carlo simulations are ~100 K lower than the MFA values for VAs and MnAs 
and ~300 K for CrAs. Note that the mean-field formula for the Curie temperature 
used in references j^Hl OH] is derived for a single-sublattice system. Therefore the 
exchange interaction between 3d and sp atoms is not treated consequently. Also the 
formula of Kiibler [52] is derived for an one-sublattice case. In general, the treatment 
of a multiple-sublattice materials needs solving a system of equations with the number 
of equations equal to the number of sublattices (section [3Jl . Below we show that for 
materials studied in this paper the one-sublattice approximation is well founded. In 
other systems, however, it leads to a large error |5"6ll57] . 

In this work we report the theoretical investigation of the magnetic exchange 
interactions for selected zb compounds. We study two Cr compounds, CrAs and 
CrSe, and four Mn compounds, MnAs, MnSi, MnGe and MnC. MnAs is isoelectronic 
to CrSe whereas other three Mn compounds are isoelectronic to CrAs. We use the 
calculated exchange parameters to estimate the Curie temperature in both the mean 
field approximation (T^f FA ) and random phase approximation (Tq PA ). A special 
attention is paid to the influence of the position of the Fermi level on the magnetic 
properties. 

The paper is organized as follows: In section we discuss the calculational 
technique. In sections [3] and 0] we present the results for the Cr and Mn based half- 
metallic compounds. Section[3]is devoted to MnC. Finally we summarize in sectional 

2. Calculational Method 

The calculations are carried out using the augmented spherical waves method (ASW) 
53 within the atomic-sphere approximation (ASA) [SI]. The exchange-correlation 
potential is chosen in the generalized gradient approximation |55j . A dense Brillouin 
zone (BZ) sampling 30 x 30 x 30 is used. The radii of all atomic spheres are chosen 
equal. 

The method of the calculation of the exchange constants has been already 
presented elsewhere [12] ■ Here, to make the paper reasonably self-contained a brief 
overview is given. We describe the interatomic exchange interactions in terms of the 
classical Heisenberg Hamiltonian 




R.R' 




RR' R R' 



(1) 



(juR^i/R') 



In equation n the indices /j, and v number different sublattices and R and R' are the 
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lattice vectors specifying the atoms within sublattices, is the unit vector pointing 
in the direction of the magnetic moment at site (/z, R). 

We employ the frozen-magnon approach 58 to calculate interatomic Heisenberg 
exchange parameters. The calculations involve few steps. In the first step, the 
exchange parameters between the atoms of a given sublattice [i are computed. 
The calculation is based on the evaluation of the energy of the frozen-magnon 
configurations defined by the following atomic polar and azimuthal angles 

0r = *. R = q-R + ^. (2) 
The constant phase is always chosen equal to zero. The magnetic moments of all 
other sublattices are kept parallel to the z axis. Within the Heisenberg model {Q the 
energy of such configuration takes the form 

(9, q) = E^ (8) + sin 2 6J^ (q) (3) 
where E^ does not depend on q and the Fourier transform J Miy (q) is defined by 

J^(q) = £j&exp(iq.R). (4) 

R 

In the case of v = \i the sum in equation @ does not include R = 0. Calculating 
E^{0 : q) for a regular q-mesh in the Brillouin zone of the crystal and performing 
back Fourier transformation one gets exchange parameters for sublattice fx. 
The determination of the exchange interactions between the atoms of two different 
sublattices fx and v is discussed in reference [ST)] . 

First, the Curie temperature is estimated within the mean-field approximation 
(MFA) for the case of a multi-sublattice material by solving the system of coupled 
equations |5TfllBTj] 

<o=3^5:w (5) 

where (s v ) is the average z component of sj^ and J(f v = Jq R - Equation [5J can be 
represented in the form of eigenvalue matrix problem 

(0 - TI)S = (6) 

where 0^ = gfj^o"^ I is a unit matrix and S is the vector of (s 1 '). The largest 
eigenvalue of matrix 8 gives the value of T^ FA [5J2] . 

Within the random phase approximation (RPA) the Curie temperature is given 
by the relation |S2] 

1 = 6/i B lr^ 1 
k B T* PA M 7V^(q)' 1 ' 

where w(q) is the spin- wave dispersion. Equation J7J) does not take into account the 
presence of several sublattices and will be used to estimate the Curie temperature 
with account for the interactions within the sublattice of the 3d atoms. As we will 
show these interactions give a leading contribution to the Curie temperature. 
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Table 1. Calculated atom-resolved and total spin moments in fiB for CrAs and 
CrSe (2 nd to 6 th column). CrAs and CrSe are half-metallic for the GaAs and CdS 
experimental lattice constants, respectively, a// means that the Fermi level is at 
the upper edge of the gap and a,ju corresponds to 1% contraction with respect 
to a//. 



Compound 


a(A) 




m As,Se 


m Voidl 


m Void2 


m Total 


CrAs - a I[GaAs] 


5.65 


3.24 


-0.31 


-0.03 


0.10 


3.00 


CrAs - a// 


5.53 


3.17 


-0.25 


-0.03 


0.11 


3.00 


CrAs - &in 


5.47 


3.10 


-0.23 


-0.03 


0.10 


2.95 


CrSe - & I[C dS] 


5.82 


3.92 


-0.10 


0.02 


0.16 


4.00 


CrSe - slji 


5.62 


3.84 


-0.04 


0.02 


0.18 


4.00 


CrSe - sun 


5.56 


3.77 


-0.04 


0.02 


0.16 


3.92 



3. Cr-compounds: CrAs and CrSe 

We begin with the discussion of the Cr compounds: pnictide CrAs and chalcogenide 
CrSe. The half-metallicity in these compounds results from the formation of bonding 
and antibonding hybrides between the t2 g states (d xy , d yz , and d xz ) of the transition 
metal atom (here Cr) and the p states of the four neighboring sp atoms (As or Se) . The 
hybridizational gap is partly filled by the e g states of the 3d atom (d z 2 and d x 2_ y 2). 
The hybridization takes place in both spin channels. The position of the hybridization 
gap is different for majority and minority electrons as a consequence of the exchange 
splitting |25| . There are in total four occupied minority states, one coming from the s 
states of the sp atom and three bonding 3d states. The minority-spin e g orbitals lay 
above the Fermi level. Thus the total spin moment M t in //b is given by the relation 
M t = [Z to t — 8) /is, where Z to t is the total number of valence electrons in the unit 
cell In the cases of CrAs (11 valence electrons) and CrSe (12 valence electrons) 
the total spin moment should be 3 and 4 respectively. 

In table^we present the atom-resolved and total spin moments. Notations Voidl 
and Void2 are used for two types of empty spheres used in our ASA calculations 
to describe correctly the charge distribution in the zinc-blende structure. For each 
system we performed calculations for three different lattice parameters. The first 
parameter (case I) corresponds to the experimental lattice constant of GaAs(CdS) for 
CrAs(CrSe). For these lattice constants both compounds are half-metals. This is seen 
in figure n where the spin-projected density of states (DOS) is presented. The total 
spin moment is exactly integer as expected for a half-metallic system. The Cr atoms 
posses a very large spin moment of 3.24 in CrAs that reaches the value of 3.92 /ib 
in CrSe. A larger moment of CrSe is a consequence of one extra valence electron per 
formula unit. Both As and Se atoms have small induced moments antiparallel to the 
Cr moments. The discussion of the spin moments of the sp atoms and their relative 
magnitude can be found in reference |25j. 

The half-metallic gap of CrSe is much larger than for CrAs. This is explained by 
the larger Cr moment in CrSe that leads to a larger exchange splitting AE X and, as 
a result, to a higher position of the minority e g bands. 

In the next step, we compress the lattice simulating a hydrostatic pressure. The 
contraction of the lattice constant results in a higher energy position of the majority 
p states with a corresponding increase of the Fermi energy, Ep. The calculations 
are performed for two decreased values of the lattice parameter. The first (case II) 
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Figure 1. Calculated spin-resolved density of states of CrAs and CrSe. The 
upper panel corresponds to the experimental lattice constant of GaAs and CdS, 
respectively, the middle panel to the aj/ case in table^ (Fermi level at the higher- 
energy edge of the half-metallic gap) and the bottom to the a.ju case in the same 
table (1% contraction of the lattice constant with respect to the a/j case). 

corresponds to the position of the Fermi level at the bottom edge of the minority-spin 
conduction band (fig. QJ. The second parameter (case III) is obtained by a further 1% 
decrease of the lattice constant. In the latter case the Fermi level lies within the spin- 
minority conduction band. Since half-metallicity is lost the total magnetic moment in 
case III becomes non-integer (table 0. 

The contraction of the lattice leads to an increase of the hybridization between 
the Cr states and the states of the sp atoms. As the result, both the Cr moment 
and the induced moment of the sp atom decrease. In case II, the changes of different 
moments compensate resulting in the same integer total moment. The decrease of the 
Cr spin moment causes a small decrease of the width of the gap. 

Next we used the frozen-magnon approach to calculate the exchange parameters 
for both systems and three values of the lattice constants (figure The exchange 
interaction between the induced moments of the sp atoms is very small and the 
corresponding parameters are not presented. 

For the semiconductor lattice parameter (left panels in figure |2J the interactions 
are characterized by the large ferromagnetic coupling of the nearest Cr atoms. In 
CrAs, this coupling is about twice the coupling in CrSe. In the zinc-blende structure 
the nearest surrounding of a Cr atom is formed by 4 sp atoms. The second-nearest 
coordination sphere contains Cr atoms. Although the exchange interaction Cr-sp atom 
is antiferromagnetic it makes the parallel directions of the Cr moments energetically 
preferable. Therefore the indirect exchange interaction of the Cr moments through 
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Figure 2. The exchange constants of CrAs (upper panel) and CrSe (lower 
panels) for three different lattice constants. The left panel corresponds to the 
experimental lattice constant of GaAs and CdS respectively for CrAs and CrSe, 
the middle one to case II and the right one to case III. 



the sp atoms stabilizes further the ferromagnetic configuration of the Cr moments (a 
similar situation was found in full-Heusler alloys |61|). The interaction is short range 
and involves only the nearest Cr atoms. 

The influence of the lattice compression is very different for the two systems. In 
CrAs, the compression leads to a small decrease of the exchange constants revealing 
the robustness of the ferromagnetism of the compound. On the other hand, in CrSe the 
effect of the compression is very strong. Already for the a// case where the system is 
still half- metallic the exchange interaction between the nearest Cr atoms is halved with 
respect to the case of the CdS lattice constant. The interaction between second- nearest 
Cr atoms strongly increases keeping antiferromagnetic character. This is reflected in 
a very small calculated Curie temperature (table When the lattice is further 
compressed and the Fermi level moves into the minority-spin conduction band, the 
ferromagnetic interaction between the nearest Cr atoms almost vanishes whereas the 
antiferromagnetic interaction between the second nearest Cr atoms increases further. 
This leads to the instability of the ferromagnetic state. 

The calculated Curie temperatures are collected in table [21 They reflect the 
properties of the exchange interactions discussed above. We report two estimations of 
the Curie temperature within the MFA, one with only the Cr-Cr interactions taken 
into account and other with both Cr-Cr and Cr-sp atom interactions considered, and 
one RPA estimation. 

In case I, for both systems the Curie temperature is much higher than the room 
temperature. For CrAs, Tq is about twice larger than in CrSe and exceeds 1000 K. 
The compression of the lattice leads in the case of CrAs to a very moderate decrease 
of the Curie temperature. In CrSe, the Tq drops strongly with decreasing lattice 
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Table 2. Calculated Curie temperatures for CrAs and CrSe. In the 3 rd and 
4 th columns the MFA-estimations of T^ FA are presented (respectively, with only 
the Cr-Cr interactions and with both the Cr-Cr and Cr-As(Se) interactions taken 
into account. In the 5 th column, the RPA-estimation of Tq is presented. Negative 
values of Tq reveal the instability of the ferromagnetic structure. 



Compound 


a(A) 


T™ FA (Cr-Cr) 


T^MFA 
1 C 


1 C 


CrAs - a.j[ G aAs] 


5.65 


1576 


1600 


1176 


CrAs - a, j j 


5.53 


1498 


1512 


1011 


CrAs - aj// 


5.47 


1376 


1387 


860 


CrSe - a I[CdS] 


5.82 


706 


709 


508 


CrSe - slji 


5.62 


109 


117 


24 


CrSe - a; a 


5.56 


-248 


-251 


-194 



parameter becoming negative in case III revealing the instability of ferromagnetism. 

Two mean-field estimations give very similar results (tableEJ) revealing the leading 
role of the Cr-Cr interactions. Therefore the neglect of the inter-sublattice Cr-sp atom 
interactions does not result in a substantial error for these compounds. 

Since the intra-sublattice interactions dominate, the use of the RPA formula JTJ) 
is justified. Overall the Tq values within RPA are by 20-30% smaller than the MFA 
ones. The difference between MFA and RPA estimations is explained by a different 
weighting of the spin- wave excitations. In the MFA all excitations are taken with 
equal weight whereas in the RPA the weight decreases with increasing energy of the 
excitation |621 163j . 

The RPA is expected to provide better estimation of the Curie temperature. In 
sectional we will discuss the relation between the MFA and RPA values of Tq in more 
detail. The only experimental information about the magnetic transition temperature 
of the systems considered is that the Curie temperature of CrAs is well above 400 K 
|12| which is in good correlation with our estimations. 

In the case of CrAs we can compare the calculated Tq values of 1600 K within 
MFA and 1176 K within RPA with the results of previous ab-initio studies. Sanyal and 
collaborators calculated the Curie temperature within MFA and found a value of 1320 
K for the GaAs experimental lattice constant and 1100 K for the InAs experimental 
lattice constant 28 . Their Monte Carlo simulations gave for the same lattice constants 
980 K and 790 K, respectively. Kiibler used a Ginzburg-Landau approach with a 
simplified form of the dynamical susceptibility and a relaxation parameter fitted to 
reproduce the experimental Curie temperature of NiMnSb. His estimation of Tq of 
CrAs gives 1041 K |S2J - Finally, Sakuma's mean-field calculation for the theoretical 
equilibrium lattice constant of 5.82 A gave Tq about 1400 K Summarizing, all 

calculations agree that the Curie temperature of CrAs should exceed or be about 1000 
K. A high Tq value that is stable with respect to the variation of the lattice constant 
makes CrAs an interesting candidate for spintronics applications. 

4. Mn-compounds: MnSi, MnGe and MnAs 

We proceed with the discussion of Mn compounds. Here we study MnSi and MnGe 
which are isoelectronic with CrAs and MnAs which is isoelectronic with CrSe. 

In tableEHwe present the atom-resolved and total magnetic moments. MnAs with 
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Compound 


a(A) 


m Mn 


m z 


m Voidl 


m Void2 


m Total 


MnAs - a J[/nP ] 


5.87 


4.16 


-0.26 


-0.01 


0.10 


4.00 


MnAs - a, u 


5.73 


4.09 


-0.20 


-0.01 


0.11 


4.00 


MnAs — a,ui 


5.68 


4.01 


-0.19 


-0.01 


0.10 


3.91 


MnGe - a J[GaAs] 


5.65 


3.43 


-0.39 


-0.08 


0.04 


3.00 


MnGe - a n 


5.61 


3.40 


-0.37 


-0.07 


0.04 


3.00 


MnGc - a,jjj 


5.55 


3.35 


-0.34 


-0.07 


0.04 


2.98 


MnSi - a, I[GaAs] 


5.65 


3.51 


-0.49 


-0.07 


0.05 


3.00 


MnSi - aj/ 


5.52 


3.42 


-0.42 


-0.06 


0.06 


3.00 


MnSi - &ju 


5.46 


3.37 


-0.39 


-0.06 


0.05 


2.97 
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Figure 3. The density of states of MnAs, MnGe and MnSi. See the caption in 
figure HI for the description of details. 



the lattice constant of InP, and MnSi and MnGe with the lattice constant of GaAs 
are half- metallic and their total spin moments following the relation presented above 
are 4 /ib for MnAs and 3 /^b for MnSi and MnGe (see tableland figure [3J ■ The Mn 
spin moment is about 3.5 /^b for MnSi and MnGe, and exceeds the 4 ^b in the case 
of MnAs reflecting the presence of one more valence electron in the latter compound. 
The moments of the sp atoms are antiferromagnetically coupled to the Mn moments. 
The value of the half-metallic gap is determined by the value of the exchange splitting 
and the energy position of the p states of the sp atom while the value of the exchange 
splitting is governed by the value of the Mn moment. As a result, the largest gap is 
obtained for MnAs followed by MnSi and the smallest gap is obtained for MnGe. 

As for the Cr compounds the contraction of the lattice increases the energy of 
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Figure 4. Variation of the exchange constants for the Mn-based compounds. 
Left panels correspond to the experimental lattice constants of semiconductors 
(see table [3J, middle panel to case II and the right panel to case III. 



the majority p states leading to an increased Fermi level. For all three compounds the 
hybridization between the p states of the sp atom and the t2 9 states of Mn increases 
and the Mn spin moment slightly decreases with the lattice contraction. The smaller 
Mn moments and, as a result, smaller exchange splittings lead to slightly smaller 
gap- widths. 

In figure 0] we present the calculated exchange constants. The exchange 
interaction between the induced moments of sp atoms is very weak and is not 
presented. There is strong similarity in the behavior of the exchange parameters for 
the Mn compounds considered here and the Cr compounds discussed in the preceding 
section. The properties of MnGe and MnSi are analogous to the properties of the 
isovalent CrAs while properties of MnAs are similar to the properties of isovalent CrSe. 
Therefore the corresponding parts of the discussion of the Cr compounds remain valid 
also for the Mn systems. 

Ferromagnetism is stabilized by both the direct Mn-Mn exchange interaction and 
by the indirect interaction through the sp atoms. The exchange interactions are short- 
range and the Mn-Mn exchange constants for MnAs are about a factor two smaller 
than for MnSi and MnGe. 

The compression has little effect on MnSi and MnGe in contrast to MnAs. Again 
there is a strong similarity between MnAs and CrSe although the dependence of 
the exchange parameters on the lattice constant in MnAs is somewhat weaker. In 
particular, the antiferromagnetic exchange interaction between second-nearest Mn 
atoms developing with contraction has lower value than in CrSe. This results 
in preserving the ferromagnetic ground state of MnAs for all lattice parameters 
considered while in CrSe the ferromagnetic state becomes unstable (table |2j . 
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Table 4. Calculated spin moments for MnAs, MnGe and MnSi. See the caption 
in table 121 for the description. 



Compound a(A) | J"MPA (Mn-Mn) T™ FA T£ PA 



MnAs 


~ a I[InP] 


5.87 


679 


701 


551 


MnAs 


- a// 


5.73 


400 


417 


275 


MnAs 


- a/// 


5.68 


201 


230 


136 


MnGe 


~ a I[GaAs] 


5.65 


1160 


1234 


716 


MnGe 


- a// 


5.61 


1180 


1241 


740 


MnGe 


- 3,111 


5.55 


1184 


1231 


695 


MnSi - 


' a-IlGaAs] 


5.65 


1255 


1402 


848 


MnSi - 


a// 


5.52 


1310 


1401 


857 


MnSi - 


&III 


5.46 


1291 


1367 


911 



The Curie temperatures of the Mn compounds are collected in tabled Again, 
as for Cr compounds, the Mn-Mn interaction is dominant and the account for the 
exchange interaction Mn-sp atom influences the value of the Curie temperature 
weakly. In agreement with the properties of the exchange constants the Curie 
temperatures of MnGe and MnSi depend weakly on the lattice parameter in contrast 
to MnAs where this dependence is strong. The Curie temperature of MnAs at highest 
contraction considered assumes the value of about ^th of the Curie temperature at 
the semiconductor lattice constant. The RPA gives values of Tq that are about 30 
% smaller than the corresponding values within the MFA. For MnSi and MnGe, also 
the RPA value of Tq remains considerably above the room temperature at all lattice 
constants studied. The stability of the high Tq value with respect to the variation 
of the lattice constant makes MnGe and MnSi promising candidates for spintronics 
applications. MnSi has an additional advantage of a larger half-metallic gap. 

Next we compare the results of our calculation of Tq with the results of previous 
calculations where available. For MnAs, we can compare our MFA-estimation of 701 
K and RPA estimation of 551 K obtained for the InP lattice constant with the Curie 
temperatures reported by Sanyal et al [2Hj ■ The MFA-estimation of Sanyal et al gave 
the value of 640 K for the InAs lattice parameter which is very close to the parameter 
used in our calculations. Their Monte Carlo simulation gave a Tq of 530 K which is 
very close to our RPA estimation. Kubler obtained the value of 671 K for MnAs with 
the InAs lattice constant. When he compressed the lattice down to the GaAs lattice 
constant of 5.65A the Tq dropped sharply to 210 K [S^j in agreement with our results. 
Finally, Sakuma's MFA-estimation of the Curie temperature of MnGe for a lattice 
constant of 5.82 A gave a value of about 1200 K which is very close to our MFA value 
of 1234 K [21 . Thus the results of different calculations are in reasonable agreement 
with one another. 

5. MnC 

The last Heusler compound studied in this paper is MnC. This system was predicted 
to be half-metallic by Qian and collaborators First we discuss the results of 

the calculation with the theoretical equilibrium lattice constant reported by Qian et 
al. The magnetic properties of MnC differ from the properties of the Mn compounds 
discussed above. MnC has 11 valence electrons: 6 of majority-spin and 5 minority- 




Figure 5. Left panel: calculated spin-resolved DOS of MnC. The upper panel 
corresponds to the lattice constant given in reference 1351 . the middle and lower 
panels to the aj/ and a,ju cases in table|S] Right panel: Atom- and spin-resolved 
DOS for the half-metallic MnC and MnSi compounds. 



spin electrons. This electron distribution results in a magnetic moment which is much 
smaller than in the compounds discussed above (table |SJ • In contrast to the other 
compounds the half-metallic gap is now in the majority-spin channel (figure [SJ). The 
majority-spin DOS has the following structure: low in energy in the interval [-0.55 Ry,- 
0.4 Ry] lie s-states. Above, in the interval [-0.4 Ry,-0.2 Ry] there are three bonding 
hybrids between the p states of C and the t2 g states of Mn. This is seen in the 
atom-projected DOS of MnC (right panel in figure |5j) where both Mn and C states 
have similar contribution in this energy interval. The group of states next in energy 
(interval [-0.2 Ry,-0.05 Ry]) does not hybridize with the lower states. These are the 
Mn e g states separated by a gap from the antibonding C p-Mn t2 g hybrids. The Fermi 
level lies within this gap in the majority-spin DOS. For the minority-spin channel, the 
position of the Fermi level belongs to the region of the e g states. These states are 
shifted by the exchange splitting to higher energies with respect to the corresponding 
majority-spin states. 

The comparison of the DOS of half-metallic MnC and MnSi (figure |SJ allows 
to reveal the reason for the move of the semiconductor-type behavior from the spin- 
minority channel in MnSi to the spin-majority channel in MnC. First, the two systems 
differ by the value of the magnetic moment that leads in the case of MnC to a smaller 
relative exchange shift of the spin-up and spin-down states and to the occupation of a 
larger number of the spin-down Mn 3d states. In particular, the minority-spin e g states 
are partly occupied in MnC whereas in MnSi these states lie above the half-metallic 
gap. Second, there is difference in the position of the e g states in the hybridizational 
gap between bonding and antibonding t2 g -p hybrids which in the case of MnSi results 
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Table 5. Calculated spin magnetic moments and Curie temperatures for MnC. 
The value of 4.20 A is the equilibrium lattice constant calculated by Qian et al. 

no 



Compound 


a(A) 


m Mn 


m z 


m Voidl 


^ Void 2 


m Total 


MnC - a J(re/ 
MnC - a/j 
MnC - a,ui 


4.20 
4.23 
4.27 


1.16 
1.16 
1.13 


-0.14 
-0.14 
-0.16 


-0.01 
-0.01 
-0.01 


-0.01 
-0.01 
-0.01 


1.00 
1.00 

1.05 


Compound 


a(A) 




T^MFA 
1 C 


(Mn-Mn) 


1 c 


jiRPA 


MnC - a J(re/ |25J) 
MnC - a// 
MnC - a/// 


4.20 
4.23 
4.27 






520 
467 
346 


526 
476 
363 


507 
458 
327 



in the gap between occupied majority e g and empty antibonding hybrid states. 

To extend the study of MnC to the states with the Fermi level lying at the top-edge 
of the half-metallic gap (case II) and above it (case III) we need in this case to expand 
the lattice (table EJ) opposite to all systems considered above where the lattice must 
be compressed. This difference comes from the move of the semiconductor channel to 
the majority-spin subsystem. The expansion of the lattice makes the energy bands 
narrower and decreases the strength of the 3d-p hybridization. At the same time the 
Mn moment increases increasing the exchange splitting between majority and minority 
states. These properties lead to increasing occupation of the majority-spin states and 
the loss of the half-metallicity (figure 0). 

The calculated exchange constants are presented in figure 0] The ferromagnetism 
is stabilized by the direct Mn-Mn interactions and the Mn-C ferromagnetic coupling. 
Since the Mn moment is much smaller than in other compounds the value of the Mn- 
Mn exchange interaction is substantially lower (figure^}. On the other hand, in MnC 
the first and second nearest Mn atoms are characterized by a similar strength of the 
ferromagnetic exchange interaction while in other systems only the nearest Mn atoms 
interact strongly. The expansion of the lattice leads to decreasing Mn-Mn exchange 
interactions. 

The calculated Curie temperatures are presented in table El The Curie 
temperature of MnC is smaller than in other Mn compounds being about 500 K 
(it is comparable with the RPA estimation of Tq of MnAs). The value of the Curie 
temperature decreases with expansion of the lattice showing again that half-metallicity 
is favorable for high Curie temperature. 

A remarkable feature of the calculation for MnC is a very small difference between 
the Tc values estimated within the MFA and the RPA. In section|2|we briefly discussed 
the reason behind the different values of the Curie temperatures estimated within MFA 
and RPA. A striking contrast between MnC and the systems considered in the previous 
sections allows gaining a deeper insight into the physics of the formation of the Curie 
temperature. We will base our analysis on the comparison of two systems: MnC and 
MnSi. Characterizing the relative difference of the MFA and RPA values of the Curie 
temperature by the relation [T<¥ FA - T^ PA ]/T^ FA we get for MnC a small value of 
2.5% compared to a large value of 47% in MnSi. 

In MFA, the Curie temperature is determined by an arithmetic average of the 
magnon energies while in RPA Tc is determined by the harmonic average of the same 
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Figure 6. Frozen-magnon energies in MnC and MnSi for [0 q] and [0 q q] 
directions in the Brillouin zone. The vertical broken line shows the Brillouin zone 
boundary for the [0 q q] direction. 



quantities. It is an arithmetic property that the MFA estimation is larger than the 
RPA one or equal to it if all averaged numbers are equal to each other. In terms of 
magnon energies, T^ FA is equal to T^ PA in the case that the magnon spectrum is 
dispersion-less: the magnon energy does not depend on the wave vector q. 

In figure El we plot the frozen-magnon energies of MnSi and MnC for two selected 
directions in the reciprocal space. We see that the MnC curves are flat in the second 
half of the q interval showing here a very weak dispersion. On the other hand, the 
MnSi curves are much closer to a simple cosinusoidal form with considerable dispersion 
in this part of the q interval. Comparing now the low-q parts of the curves we notice 
that MnC curves lie higher than the MnSi curves. Therefore, the relative contribution 
of the low-energy magnons to the RPA value of the Curie temperature is less important 
in MnC compared to MnSi. 

This combination of features of the wave-vector dependencies of the frozen- 
magnon energies is responsible for a much smaller difference between the RPA and 
MFA estimations of the Curie temperature in the two cases. Since the form of 
the magnon dispersion is reflected on the properties of the inter-atomic exchange 
interactions it is instructive to compare the patterns of the exchange interactions 
of the two systems. Taking into account that only Mn-Mn interactions contribute 
importantly in the Curie temperature we will focus on the comparison of these 
interactions. 

The analysis of the Mn-Mn exchange interactions in MnSi and MnC (figure 
shows a strong difference between the two systems. In MnSi, only the first-neighbor 
exchange interaction plays important role. All other interactions are much smaller. 
On the other hand, in the case of MnC the second-neighbor interaction is comparable 
with the first-neighbor one. The further few interactions though sizable play less 
important role. The presence of two first coordination spheres with similar values of 
the exchange interactions, on the one hand, makes the exchange interaction in MnC 
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more long ranged and, on the other hand, it effectively increases the coordination 
number of the atoms with leading exchange interaction. Either of these properties 
leads to a better accuracy of the MFA estimation. 

6. Summary and Conclusions 

We have presented a detailed study of the stability of the ferromagnetism in six 
transition-metal compounds crystallizing in the zinc-blende structure. We calculated 
from first-principles the electronic structure and employed the frozen-magnon 
approximation to evaluate the exchange parameters and the Curie temperature. In the 
calculation of the Curie temperature we used both the mean-field and random-phase 
approximations . 

Considering the CrAs and CrSe compounds we found in the case of CrAs a very 
high Curie temperature exceeding 1000 K. The ferromagnetism of CrAs is stable with 
respect to the compression of the lattice. On the other hand CrSe possesses about 
two times lower Curie temperature which decreases strongly with lattice compression. 

The properties of MnSi and MnGe are similar to the properties of the isovalent 
CrSe while the properties of MnAs are similar to the properties of isovalent CrSe. 

MnC presents an interesting case that differs from all other systems studied. This 
system is also half-metallic with the half-metallic gap in the majority-spin channel in 
contrast to other systems where the half-metallic gap is in the minority-spin subsystem. 
The Curie temperature of MnC was found to be about 500 K. 

In all cases considered we found that half-mctallicity is favorable for 
ferromagnetism and leads to an increased Curie temperature. 

We conclude that two compounds, MnSi and CrAs, possess properties which 
make them promising candidates for spintronics applications: large half-metallic gap 
and high Curie temperature. The stability of the Curie temperature in these systems 
allows to expect that being grown on the substrate of a binary semiconductor they 
will preserve the desired magnetic properties. 

Acknowledgments 

The financial support of Bundesministerium fiir Bildung und Forschung is 
acknowledged. 

References 

[1] Prinz G A 1998 Science 282 1660 

[2] de Boeck J, van Roy W, Das J, Motsnyi V, Liu Z, Lagae L, Boeve H, Dessein K and Borghs G 

2002 Semicond. Sci. Technol. 17 342; de Boeck J, van Roy W, Motsnyi V, Liu Z, Dessein K 

and Borghs G 2002 Thin Solid Films 412 3 
[3] de Groot R A, Mueller F M, van Engcn P G and Buschow K H J 1983 Phys. Rev. Lett. 50 2024 
[4] Soulen Jr R J, Bycrs J M, Osofsky M S, Nadgorny B, Ambrose T, Cheng S F, Broussard P R, 

Tanaka C T, Nowak J, Moodera J S, Barry A and Cocy J M D 1998 Science 282 85 
[5] Schwarz K J. Phys. F: Met. Phys. 1986 16 L211; Korotin M A, Anisimov V I, Khomskii D I 

and Sawatzky G A 1998 Phys. Rev. Lett. 80 4305 
[6] Yanase A and Siratori K 1984 J. Phys. Soc. Jpn 53 312 

[7] Galanakis I, Dederichs P H and Papanikolaou N 2002 Phys. Rev. B 66 134428 
[8] Galanakis I, Dederichs P H and Papanikolaou N 2002 Phys. Rev. B 66 174429 
[9] Galanakis I 2004 J. Phys.: Condens. Matter 16 3089 
[10] Akai H 1998 Phys. Rev. Lett. 81 3002 

[11] Sandratskii L M and Bruno P 2002 Phys. Rev. B 66 134435 



Stability of ferromagnetism in the half-metallic pnictides 



16 



[12] Akinaga H, Manago T and Shirai M 2000 Japan. J. Appl. Phys. 39 L1118 
[13] Shirai M 2001 Physica E 10 143 

[14] Shirai M 2003 J. Appl. Phys. 93 6844; Shirai M 2004 J. Phys.: Condens. Matter 16 S5525 
[15] Mizuguchi M, Akinaga H, Manago T, Ono K, Oshima M and Shirai M 2002 J. Magn. Magn. 
Mater. 239 269 

[16] Zhao J H, Matsukura F, Takamura K, Abe E, Chiba D and Ohno H 2001 Appl. Phys. Lett. 79 
2776 

[17] Ono K, Okabayashi J, Mizuguchi M, Oshima M, Fujimori A and Akinaga H 2002 J. Appl. Phys. 
91 8088 

[18] Mizuguchi M, Akinaga H, Manago T, Ono K, Oshima M, Shirai M, Yuri M, Lin H J, Hsieh H 

H and Chen C T 2002 J. Appl. Phys. 91 7917 
[19] Landolt-Bornstein, New Series, Group III, Vol. 27, Pt.a edited by K Adachi and S Ogawa 

(Springer- Verlag, Heidelberg) 1988 
[20] Wei S H and Zunger A 1987 Phys. Rev. B 35 2340 
[21] Sanvito S and Hill N A 2000 Phys. Rev. B 62 15553 

[22] Continenza A, Picozzi S, Ceng W T and Freeman A J 2001 Phys. Rev. B 64 085204 
[23] Zhao Y -J, Ceng W T, Freeman A J and Delley B 2002 Phys. Rev. B 65 113202 
[24] Xie W -H, Xu Y-Q, Liu B G and Pettifor D G 2003 Phys. Rev. Lett. 91 037204 
[25] Galanakis I and Mavropoulos Ph 2003 Phys. Rev. B 67 104417 
[26] Liu B G 2003 Phys. Rev. B 67 172411 

[27] Pask J E, Yang L H, Fong C Y, Pickett W E and Dag S 2003 Phys. Rev. B 67 224420 

[28] Sanyal B, Bergqvist L and Eriksson O 2003 Phys. Rev. B 68 054417 

[29] Xu Y G, Liu B G and Pettifor D G 2002 Phys. Rev. B 66 184435 

[30] Xie W -H and Liu B G 2003 J. Phys.: Condens. Matter 15 5085 

[31] Yao K L, Gao G Y, Liu Z L and Zhu L 2005 Solid State Commun. 133 301 

[32] Zhang M, Hu H, Liu G, Y. Cui, Liu Z, Wang J, Wu G, Zhang X, Yan L, Liu H, Meng F, Qu 
J and Li Y 2003 J. Phys.: Condens. Matter 15 5017; Zhang M, Hu H, Liu G, Liu Z, Cui Y 
and Wu G 2004 J. Low Temp. Phys. 135 267; Zhang M, Hu H, Liu Z, Liu G, Cui Y and Wu 
G 2004 J. Magn. Magn. Mater. 270 32 

[33] Xie W -H, Liu B G and Pettifor D G 2003 Phys. Rev. B 68 134407 

[34] Zhang M, Briick E, de Boer F R, Liu G D, Hu H N, Liu Z H, Cui Y T and Wu G C 2004 J. 

Mater. Res. 19 2738 
[35] Qian M C, Fong C Y and Yang L H 2004 Phys. Rev. B 70 052404 
[36] Nakao M 2004 Phys. Rev. B 69 214429 

[37] Qian M C, Fong C Y, Pickett W E and Wang H Y 2004 J. Appl. Phys. 95 7459 
[38] Shi L J and Liu B G 2005 J. Phys.: Condens. Matter 17 1209 

[39] Yamana K, Geshi M, Tsukamoto H, Uchida I, Shirai M, Kusakabc K and Suzuki N 2004 J. 

Phys.: Condens. Matter 16 S5815 

[40] Zhao Y - J and Zunger A 2005 Phys. Rev. B in press; preprint |arXiv:conqTm at /0408707 
[41] Zheng J -C and Davenport J W 2004 Phys. Rev. B 69 144415. 
[42] Galanakis, I 2002 Phys. Rev. B 66 012406 

[43] Fong C Y, Qian M C, Pask J E, Yang L H and Dag S 2004 Appl. Phys. Lett. 84 239; Fong C Y 

and Qian MC 2004 J. Phys.: Condens. Matter 16 S5669. 
[44] Mavropoulos Ph and Galanakis I 2004 J. Phys.: Condens. Matter 16 4261 
[45] Nagao K, Shirai M and Miura Y 2004 J. Appl. Phys. 95 6518 
[46] Geshi M, Kusakabe K and Suzuki N 2004 J. Phys.: Condens. Matter 16 S5701 
[47] Qian M C, Fong C Y, Pickett W E, Pask J E, Yang L H and Dag S 2005 Phys. Rev. B 71 012414 
[48] Bengone O, Eriksson O, Fransson J, Turek I, Judrnovsky J and Drchal V 2004 Phys. Rev. B 70 

035302 

[49] Mavropoulos Ph, Sato K, Zeller R, Dederichs P H, Popescu V and Ebert H 2004 Phys. Rev. B 
69 054424 

[50] Chioncel L, Katsnelson M I, de Wijs G A, de Groot R A and Lichtenstein A I 2005 Phys. Rev. 
B 71 085111 

[51] Sakuma A 2002 J. Phys. Soc. Japan 71 2534 
[52] Kiibler J 2003 Phys. Rev. B 67 220403 

[53] Williams A R, Kiibler J and Gclatt C D 1979 Phys. Rev. B 19 6094 

[54] Andersen O K 1975 Phys. Rev. B 12 3060 

[55] Perdew J P and Wang Y 1992 Phys. Rev. B 45 13244 

[56] §a§ioglu E, Sandratskii L M and Bruno P 2004 Phys. Rev. B 70 024427 

[57] §as,ioglu E, Sandratskii L M and Bruno P 2005 J. Phys.: Condens. Matter 17 995 

[58] Rosengaard N M and Johansson B 1997 Phys. Rev. B 55 14975; Halilov S V, Eschrig H, Perlov 



Stability of ferromagnetism in the half-metallic pnictides 



17 



A Ya and Oppeneer P M 1998 Phys. Rev. B 58 293; Sandratskii L M and Bruno P 2003 Phys. 
Rev. B 67 214402 

[59] Anderson P W, 1963, Theory of magnetic exchange interactions: Exchange in insulators and 
semiconductors, in Solid State Physics, edited by F. Seitz and Turnbull (Academic Press, 
New York), Vol. 14 pp. 99-214. 
[60] §a§ioglu E, Sandratskii L M and Bruno P 2005 J. Magn. Magn. Mater. 290-291 385. 
[61] Kiibler J, Williams A R and Sommers C B 1983 Phys. Rev. B 28 1745 
[62] Pajda M, Kudrnovsky J, Turek I, Drchal V and Bruno P 2001 Phys. Rev. B 64 174402 
[63] Bouzerar G, Kudrnovsky J, Bergqvist L, and Bruno P, 2003 Phys. Rev. B 68, 081203 



